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ABSTRACT 

Optical spectra were obtained of the optical counterpart of the high latitude (6 ~ 62°) soft 
X-ray transient XTE J1118+480 near its quiescent state (i? ~ 18.3) with the new 6.5 m MMT 
and the 4.2 m WHT. The spectrum exhibits broad, double-peaked, emission lines of hydrogen 
(FWHM ~ 2400 km s^^) arising from an accretion disk superposed with absorption lines of a late- 
type secondary star. Cross-correlation of the 27 individual spectra with late-type stellar template 
spectra reveals a sinusoidal variation in radial velocity with amplitude K — 701 ± 10 km s^^ and 
orbital period P = 0.169930 ± 0.000004 d. The mass function, 6.1 ± 0.3 Mq, is a firm lower limit 
on the mass of the compact object and strongly implies that it is a black hole. We estimate the 
spectral type of the secondary to be K7V-M0V and that it contributes 28 ± 2% of the light in 
the 5800-6400 A region on 2000 November 20 increasing to 36 ± 2% by 2001 January 4 as the 
disk fades. Photometric observations (i?-band) with the lAC 0.8 m telescope reveal ellipsoidal 
light variations of full amplitude 0.2 mag. Modeling of the light curve gives a large mass ratio 
~ 20) and a high orbital inclination {i = 81° ±2°). Our combined fits yield a mass of the 
black hole in the range Mi = 6.0— 7.7 Mq (90% confidence) for plausible secondary star masses 
of M2 = 0.09 — 0.5 Mq. The photometric period measured during the outburst is 0.5% longer 
than our orbital period and probably reflects superhump modulations as observed in some other 
soft X-ray transients. The estimated distance is d = 1.9 ± 0.4 kpc corresponding to a height of 
1.7 ± 0.4 kpc above the Galactic plane. The spectroscopic, photometric, and dynamical results 
indicate that XTE Jll 18-1-480 is the first firmly identified black hole X-ray system in the Galactic 
halo. 
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1. Introduction 

X-ray novae or soft X-ray transients comprise 
a subset of low-mass X-ray binaries (LMXBs) 
which consist of a late-type secondary star and 
a neutron star or a black hole exhibiting bright 
optical and X-ray outbursts which are recurrent 
on time scales of decades (Tanaka & Lewin 1995; 
Tanaka & Shibazaki 1996). During their out- 
bursts, they resemble persistent LMXBs in which 
the light of the secondary star is overwhelmed by 
a luminous accretion disk surrounding the com- 
pact object. After a year or less in some objects, 
the system returns to quiescence. The secondary 
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star now contributes a much larger fraction of the 
total light and its atmospheric absorption lines 
become visible in optical spectra. Thus quies- 
cent X-ray novae provide the ideal opportunity 
to study the nature and dynamical properties of 
the binary system (see e.g. Charles 1998). These 
studies have demonstrated that the mass of the 
compact object in ten X-ray novae (Charles 1998; 
McClintock 1998; Nova Vel 1993: Filippenko et 
al. 1999; V4641 Sgr: Orosz et al. 2001) ex- 
ceeds the theoretical maximum mass of a neutron 
star (Rhoades & Ruffini 1974) and thus must evi- 
dently be a black hole. 

A previously unknown X-ray transient, XTE 
Jl 118-1-480, was discovered by the RXTE ASM in- 
strument on 2000 March 29 (Remillard et al. 2000) 
with an average 2-12 keV X-ray intensity of 
only 39 mCrab, but with a spectral signature 
similar to other accreting Galactic black holes. 
An optical counterpart was identified on March 
30 (Uemura et al. 2000) with V = 12.9 and 
confirmed spectroscopically (Garcia et al. 2000). 
Photometry during the outburst (Patterson 2000; 
Uemura et al. 2000) revealed a weak modulation 
with a period of~ 0.1708 d and an amplitude 
of 0.055 mag reflecting the probable orbital pe- 
riod of the system. The shape of the light curve 
and its temporal evolution resembled those of su- 
perhumps (Uemura et al. 2000) observed during 
superoutbursts of short-period cataclysmic vari- 
ables and outbursts of some other soft X-ray tran- 
sients (O'Donoghue & Charles 1996), and which 
correspond to a period only fractionally (<l-2%) 
longer than the orbital period. 

By 2000 August 31, photometry reported to 
VSNET indicated that Jl 118+480 had faded to 
near its quiescent level of y ~ 18.8. Its high 
Galactic latitude (b = +62°), implying a low value 
of interstellar absorption (Av= 0.04 mag; Hynes 
et al. 2000), combined with its relatively bright 
quiescent magnitude makes J1118+480 ideal for 
detailed study in quiescence, especially at wave- 
lengths not easily observable in other low-latitude 
transients. In this Letter, we present the results 
of optical spectroscopy and photometry obtained 
in quiescence beginning in late 2000 November 
and continuing through 2001 January. A prelim- 
inary announcement of our spectroscopic results 
(Wagner et al. 2000) as well as those obtained by 
a second group (McClintock et al. 2000) was re- 



ported in the lAU Circulars. 

2. Observations 

We obtained 20 spectra of J1118+480 with the 
newly upgraded 6.5 m MMT equipped with the 
blue channel CCD spectrograph on the nights of 

2000 November 20 (5 spectra), 30 (9 spectra), and 

2001 January 4 (6 spectra). The spectra cover 
the range AA4200-7500 at a spectral resolution of 
3.5 A and a dispersion of 1.1 A pixel" ^. The see- 
ing was typically 1" and a 1" wide slit was em- 
ployed. Each 1440 s exposure was bracketed by 
a HeNeAr lamp spectrum which led to a wave- 
length calibration accurate to 5-7 km s^^ rms. 
Spectra of the late-type dwarf stars BD+63°137 
(K7), Gliese 239 (MO), Gliese 96 (MO. 5), Gliese 
154 (MO), and Gliese 388 (AI3.5) were obtained 
using the same instrumental configuration so as 
to accurately gauge the spectral type of the sec- 
ondary star and to provide templates for the radial 
velocity analysis. We note that the first two stars, 
BD+63°137 and Gliese 239, are classified as K5V 
and K7 respectively in Simbad. However, based 
on the depth of their TiO bands and by compar- 
ison with very securely classified templates (such 
as 61 Cyg A & B), we support classifications of 
K7V and MOV respectively. 

We also obtained 7 spectra of J1118+480 with 
the ISIS red channel of the 4.2 m William Herschel 

Telescope on La Palma on the night of 2001 Jan- 
uary 11. We employed the R316R grating which 
yields 1.47 A pixer^ in the range AA5820-7320. 
The slit width was 1'.'5 which resulted in a spec- 
tral resolution of 4.5 A. Exposure times were 1200 
s and the wavelength calibration provided by in- 
ternal lamps was checked with respect to night- 
sky emission lines to be within 10 km s~^. Spec- 
tra of the dwarf stars HR 5265 (K3), HR 8085 
(61 Cyg B; K5), HR 8086 (61 Cyg A; K7), and 
Gliese 361 (Ml. 5) were both employed as radial 
velocity templates and for spectral classification 
of the secondary. Accurate absolute photometric 
calibration was not attempted for either the MMT 
or WHT spectra due to slit losses. 

Time-resolved differential CCD photometry 
was obtained of J1118+480 in the i?-band (Ac 
= 6500 A; FWHM = 1500 A) on the nights of 
2000 December 14. 28, and 2001 January 9 (UT) 
with the 0.80 m IAC80 telescope at the Observa- 
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torio del Tcidc (Tcncrifc) . We used the Thompson 
1024 X 1024 pixel CCD camera with 2x2 pixel 
binning giving an image scale of 0'/43 pixel" ^. A 
total of 68 10-mimitc exposures were obtained in 
~ 2" seeing and good transparency, together with 
bias and flatfield frames to facilitate the standard 
data reduction procedures. We applied profile fit- 
ting photometry to Jl 118+480 and several nearby 
comparison stars using IRAF. We also selected 
comparison stars which were checked for variabil- 
ity during the night and calibrated as field refer- 
ence stars using Landolt standard stars to place 
our photometry on an absolute scale. We esti- 
mate that our relative photometric accuracy is 
0.05 mag. Finally, we measured the position of 
J1118+480 on a registered average of the frames 
obtained on December 28. Based on a grid of 
22 USNO A2.0 stars, we find that the position 
of J1118+480 is R.A. = ll'*18"10!84, DecL = 
-|-48°02'12'.'9 (equinox 2000.0, accuracy iO'.'a). 

3. Results 

3.1. Spectroscopy 

The individual spectra of J1118-I-480 exhibit 
broad Balmer emission lines arising from an ac- 
cretion disk superposed with absorption lines or 
bands of Mg6 A5175, Na D, and of TiO. Radial 
velocities were extracted by cross-correlating the 
individual normalized spectra over the range 5800- 
6400 A of J1118-h480 with our grid of templates 
using routines in both IRAF and locally developed 
software packages. Phasing of the radial veloci- 
ties on the photometric period (0.17 d) yielded a 
nearly sinusoidal variation in radial velocity indi- 
cating that the photometric period was close to 
the orbital period. 

To obtain an accurate measurement of the or- 
bital period and determine the best-fitting spec- 
tral type of the secondary, successive circular 
orbits with periods ranging from 0.169000 to 
0.171000 d were fit to the radial velocity data 
defined by the grid of templates. We found that 
our best cross-correlations and circular orbit solu- 
tions corresponded to spectral types of K5V-M0V 
and an orbital period of 0.169930 ± 0.000004 d 
(;^2 ^ 4 5. 23 dof), i.e. 0.5% shorter than the out- 
burst photometric period of Uemura et al. (2000). 
We have scaled the velocity errors by a factor of 
\fAj) in order to give = l-O- 



The radial velocity curve of the secondary star 
is shown in Figure 1 (upper panel). A sine fit 
to the data gives: P = 0.169930 ± 0.000004 d; 
7 = -15±10km s"^ K = 701 ±10 km s'^To = 
HJD 2,451,868.8916 ± 0.0004 (2000 November 
20.3916 UT), where Tq corresponds to absolute 
phase 0.0 or inferior conjimction of the secondary 
star. Orbital smearing will systematically under- 
estimate the true velocity amplitude. For these 
data at the peak of the radial velocity curve, we 
estimate this effect to be no more than 1.6% or 
about 11 km s~^. The mass function is 

where G is the Gravitational Constant. The mass 
function implies a lower limit to M\ of M\ = 
6.1 Mq. This greatly exceeds the maximum mass 
of a neutron star (Rhoades & RufRni 1974) and 
implies that the compact object is probably a 
black hole. 

In Figure 2, we show the spectrum of Jl 118-1-480 

on 2000 November 30, formed by averaging the 
individual Doppler-corrected spectra (total inte- 
gration time 3.6 hr) to the rest frame of the sec- 
ondary star, together with the spectra of a M0.5V 
(Gliese 96) and K7V (BD-h63 137) star for com- 
parison, bracketing the likely spectral type of the 
J1118-t-480 secondary. The absorption features of 
the secondary star are now sharp and readily ap- 
parent, but this process smears out the emission 
lines arising from the accretion disk surround- 
ing the compact object. However, comparison of 
the spectrum with our templates indicates that 
the spectrum of Jl 118-1-480, apart from the emis- 
sion lines, is too blue. Also, the absorption lines 
are much weaker relative to an uncontaminated 
K7V star, thereby suggesting a substantial con- 
tribution from the quiescent accretion disk. This 
contribution, particularly the velocity information 
in the broad and complex Ha emission line, will 
be discussed in a forthcoming companion paper 
(Zurita et al. 2001). 

To estimate the fraction of the observed flux 
arising from the secondary star and disk sep- 
arately, we subtracted fractions of our grid of 
templates from the Doppler-corrected spectra of 
J1118+480 for each night and performed a 
minimization of the residuals in the range 5800- 
6400 A. This region overlaps both the MMT and 
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WHT spectra and is rich in metallic absorption 
lines. The templates were broadened by convolu- 
tion with Gaussian passbands to account for small 
differences in instrumental resolution between the 
MMT and WHT and the effects of orbital smear- 
ing due to the length of our exposures. We find 
that the best-fitting templates are K7V-M0V and 
that the secondary contributes ~ 32% of the total 
light in this spectral region. Our results are con- 
sistent with those obtained by McClintock et al. 
(2001) close to our earlier observations. Specif- 
ically, we find that for 2000 November 20, 30, 
and 2001 January 4 that the secondary contributes 
28 ± 2%, 28 ± 2%, 36 ± 2% respectively, so that it 
appears that the disk luminosity is still declining 
and the system may not yet be fully quiescent. 

3.2. Photometry 

To interpret the optical light curves, we used a 

model that includes a Roche lobe-filling secondary 
star, a concave accretion disk, and mutual eclipses 
of the disk and the secondary star. This X-ray 
binary model and the fitting procedure are fully 
described in Shahbaz et al. (2001). The model 
parameters are the binary inclination i, the mass 
ratio q {=Mi/M2), the mean temperature T and 
gravity log^ of the secondary, the gravity dark- 
ening exponent (3, the accretion disk radius Rdisk 
[defined as the fraction of the distance to the inner 
Lagrangian point (Rli)], the flare angle of the ac- 
cretion disk (9), the temperature at the outer disk 
edge T(iisk , and the exponent on the power-law ra- 
dial temperature distribution rj. We fixed Ay at 
0.04 mag (Hynes et al. 2000). 

From our optical spectroscopy, the secondary 
star has an observed spectral type of K7V-M0V, so 
we fixed T at 4250 K and log 5 at 5.0, appropriate 
for K7V. Since the late-type star is convective, we 
take P = 0.08 (Lucy 1967). The accretion disk is 
assumed to be optically thick and geometrically 
thin, and so we set the flare angle of the disk to 
be 0.04 radians and rj to be -0.75 (appropriate for 
a steady-state disk; Pringle 1981). We can use 
the value for the mass function and the observed 
spectral type to place an upper limit to the binary 
mass ratio. The mass function gives a lower limit 
to Ml (> 6.1 Mq) and the mass for the main 
sequence spectral type of the secondary gives an 
upper limit for M2 (< 0.52 Mq). Thus we obtain 
a lower limit for the mass ratio of q > 12 and 



so we assume g — 20 in our model (the results 
are largely insensitive to this value). Note that 
such high mass ratios are also supported by our 
small "superhump period" excess in the model of 
Mineshige et al. (1992). 

The phase-averaged i?-band light curve of 
J1118-I-40 is shown in the lower panel of Figure 
1 together with our best-fitting model for the pa- 
rameters i, Rdisk, Tdisk, a phase shift,the distance 
d, and normalization. The data were folded using 
the spectroscopic period and epoch given above 
(i.e. phase is defined as superior conjunction 
of the compact object; note that this is not how 
McClintock et al. (2001) have defined their phase 
convention). We find a best fit at i = 81 ± 2 
degrees, Rdisk = 0.8-Rli, Tdisk = 3906 K, and 
d = 1.8 ± 0.3 kpc {2a). A one-dimensional grid 
search was performed to determine the uncertain- 
ties of the fitted parameters. The model predicts a 
veiling of 76% which is consistent with the values 
derived from our optical spectroscopy. It should 
be noted that with such a high inclination one 
might normally expect to see eclipses of the sec- 
ondary and accretion disk. However, given the ex- 
treme mass ratio, one only expects shallow eclipses 
of the secondary and disk (Figure 1). 

In addition, we find a small phase shift (0.02 ± 
0.04; -^4.6 min) between the photometric light 
curve and the spectroscopically defined phase. 
Such an effect (at the level of 12 ± 2 min) is seen by 
McClintock et al. (2001) who, as we do, discount 
the possibility of an instrumental origin. It should 
be noted that McClintock et al. do not formally fit 
their data with an ellipsoidal model, but show one 
that "matches" their data, albeit with a very obvi- 
ous phase offset of 12 min. Their model combines 
an ellipsoidal modulation with a constant contri- 
bution (at the 66% level) from an accretion disk. 
We believe that this phase shift is actually due 
to a temporally varying, asymmetrically emitting 
disk component. In a companion paper (Shahbaz 
et al. 2001), we model the quiescent light curves 
of the neutron star SXT XTE .12123-058 (it has 
a 6 hr orbital period and high orbital inclination) 
which are remarkably similar in shape to that of 
McClintock et al. (2001) by employing such a disk 
component. 
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4. Discussion 

The mass of the secondary star can be con- 
strained by making assumptions regarding its evo- 
lutionary state. An upper Hmit to M2 can be ob- 
tained by using the mass of a main sequence star of 
the same spectral type (M2(MS')), while the sec- 
ondary would not have left the main sequence if its 
mass were less than the Schonbcrg-Chandrasekhar 
limiting value of 0.17M2(M5'). Thus from the ob- 
served spectral type of K7V, M2 must lie in the 
range 0.09-0.50 Mq. Combining this constraint 
with our measured values for f{M) and i, we show 
in Figure 3 the inferred range for Mi which is 
6.0-7.7 Mq {2a). The uncertainties were deter- 
mined by Monte Carlo simulation, in which we 
draw Gaussian-distributed random values for the 
observed quantities, with a mean and variance the 
same as the observed values. 

Over the past 10 years dynamical studies have 
established that several X-ray novae contain com- 
pact objects with masses significantly larger than 
the maximum mass of a normal neutron star (~ 
3.2 Mq; Rhoades & Ruffini 1974) such as V404 
Cyg (> 6.1 M©; Casares and Charles 1994); QZ 
Vul 1988 (> 5.0 Mq- Casares, Charles, & Marsh 
1995; Filippenko, Matheson, & Earth 1995; Har- 
laftis, Hornc, & FiHppenko 1996); and V2107 Oph 
1977 (> 4.7 Mq; Remillard et al. 1996; Filip- 
penko et al. 1997; Harlaftis et al. 1997) implying 
that they are black holes. The large mass function 
and well-determined mass of the compact object 
in XTE J1118-F480 derived here indicates that it 
almost certainly contains a black hole rather than 
a neutron star. Recently Bailyn et al. (1998) have 
noted that the masses of black holes in LMXBs 
are strongly peaked at ^7M0, and our mass de- 
termination is consistent with this result. 

We can estimate the distance to Jl 118-1-480 by 
constraining the secondary's size from our pho- 
tometric modeling by using the Eaton & Foe 
(1984) redetermination of the Barnes-Evans rela- 
tion, which relates y-band flux at the star and 
color for stars with measured angular diameters. 
Prom our light curve the mean R is 18.17, of which 
76% is from the disk according to our model fits. 
This implies that the secondary has R = 19.72, 
and if its spectral type is K7V {V — R = 1.15; 
from our spectroscopy), then Vq = 20.83 (assum- 
ing Av= 0.04). From Vq and the V — R color. 



we then estimate that the angular diameter of the 
secondary is 1.94 x 10~^ mas (Eaton & Poe 1984). 
Prom the Eggleton (1983) approximation of the 
volume radius of the Roche lobe of the secondary 
star, we find that the size of the Roche-lobe is 
i?2 = 0.17o for g = 20 and where a is the separa- 
tion of the center of mass of the two components. 
From our total mass estimate of 7 Mq and Ke- 
pler's third law, we find that a = 2.37Rq, and 
thus d = 1.9 ± 0.4 kpc. The error in d is domi- 
nated by the uncertainty in the secondary spectral 
type. We find that for b = +62°, z = 1.7±0.4 kpc, 
where z is the vertical distance above the Galactic 
plane. 

Thus, XTE J1118-F480 is the first firmly identi- 
fied black hole X-ray binary system in the Galac- 
tic halo. This fact is quite remarkable compared 
with the distribution of the black-hole LMXBs 
(White & van Paradijs 1995) where the mean z is 
around 400 pc (and that is dominated by a single 
object, H1705-25, at ~950 pc). This is less than 
half that of the neutron star systems, which is ex- 
plained by the additional "kick" velocity received 
in the formation of the neutron star. To place 
Jl 118-1-480 at its current location in this scenario 
would have required an enormous velocity out of 
the plane, making its almost zero systemic veloc- 
ity very unlikely. Furthermore, we find no evi- 
dence of a large proper motion of J1118-(-480 (< 10 
mas yr~^;2iT) since its position measured during 
the outburst (Uemura et al. 2000; Masi 2000) and 
our position near quiescence (see §2) as compared 
with its position measured on an archival POSS-I 
plate obtained ~ 47 yrs ago (from the USNO A2.0 
and APM Palomar Schmidt sky catalog; Lewis & 
Irwin 1996) all agree to within the astrometric er- 
rors. This object therefore presents interesting 
challenges and constraints on the formation and 
evolution of SXTs. 

It is also interesting to consider whether there 
might be an additional population of such high z 
BH LMXBs in the Galaxy. J1118-h480 is one of 
the closest, yet is also one of the (intrinsically) 
X-ray weakest of the LMXB X-ray transients and 
has a low Lx/Lopt ratio. Hence, further exam- 
ples would almost certainly have gone undetected 
due to their X-ray faintness and high latitude (the 
most sensitive X-ray monitoring is undertaken in 
the plane and around the Galactic center). How- 
ever, Hynes et al. (2000) discussed the possibility 
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that the intrinsic weakness of the X-rays relative 
to the optical might be due to the high inclination, 
making this system akin to the Galactic ADC (ac- 
cretion disk corona) sources, and hence other high 
latitude systems would be expected to be much 
brighter X-ray sources. While this appears to have 
been borne out by our inference of the (extremely) 
high inclination of 81°, there are several difScul- 
ties with this argument. There was an absence 
of any X-ray modulation in outburst, although 
the secondary is so small that it might always be 
completely shadowed by the disk and hence un- 
able to produce any significant modulation. (Note 
that this interpretation is testable via higher time 
resolution photometry as the eclipse should be 
visible; see Figure 1 at phase 0.0). More seri- 
ously, the models considered by Esin et al. (2001) 
and fit to the exceptionally wide wavelength spec- 
tra of J1118+480 obtained during outburst have 
ruled out the ADC model since the source was 
in the low/hard state. This is further supported 
by the direct X-ray observation of ^O.lHz QPOs 
(Wood et al. 2000). 
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Orbital Phase (P = 0.169930 d) 



Fig. 1. — Upper: Folded radial velocities of XTE 

Jl 118+480 and the best- fitting sinusoid. Open 
circles indicate the data obtained on 2000 Novem- 
ber 20 (MMT), filled circles on 2000 November 30 
(MMT), filled squares on 2001 January 4 (MMT), 
and open squares on January 12 (WHT) respec- 
tively. Lower: The phase folded i?-band light 
curve of XTE J1118+480. The solid hnc shows 
the best-fit light curve solution. The sharp dip at 
phase 1.0 might be due to a grazing eclipse of the 
accretion disk by the secondary star. Two orbital 
cycles are shown for clarity. 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 
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Fig. 2. — Upper: Doppler-corrected rest-frame 
spectrum of XTE Jl 118+480 obtained on 2000 
November 30. Middle: M0.5V spectrum (Gliese 
96). Lower: K7V spectrum (BD+63°137). We 
estimate that the secondary is K7V-M0V in 
J1118+480. The spectra have been offset for clar- 
ity. The prominent spectral features are indicated. 




Fig. 3. — The mass of the compact object as a 
function of the secondary star mass. The two solid 
lines represent the 90% confidence limits for the 
probable mass range (the central dashed line being 
the best fit value). The two dashed vertical lines 
are upper and lower limits to the secondary star 
mass from other considerations (see text). 



8 



